Effects of decoherence and errors on Bell-inequality violation 
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We study optimal conditions for violation of the Clauser-Horne-Shimony-Holt form of the Bell 
inequality in the presence of decoherence and measurement errors. We obtain all detector config- 
urations providing the maximal Bell inequality violation for a general (pure or mixed) state. We 
consider local decoherence which includes energy relaxation at the zero temperature and arbitrary 
dephasing. Conditions for the maximal Bell-inequality violation in the presence of decoherence are 
analyzed both analytically and numerically for the general case and for a number of important 
special cases. Combined effects of measurement errors and decoherence are also discussed. 
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I. INTRODUCTION 

Entanglement, i.e., quantum nonlocal correlations be- 
tween physical systems, is not only a basic feature of 
quantum behavior [l|, [2j , but also an important resource 
for quantum computation and quantum information [3j. 
Decoherence, i.e., loss of coherence of states of quan- 
tum systems due to the interaction with the environ- 
ment, is one of the major stumbling blocks for quan- 
tum computation [3[. Therefore recently there has been 
a surge of interest in effects of decoherence on entangle- 
ment H i, S 0, 1, S M EH, E3, 0, • 

One of the most striking manifestations of the non- 
classical nature of entanglement is violation of the Bell 
inequality (BI) [15]. Effects of decoherence on the BI 
violation have attracted a significant interest recently 
0, S, B M, EH, El ES EE d . Decoherence transforms a 
pure entangled state into a mixed state, decreasing thus 
entanglement and the Bell inequality violation. While 
any pure (completely coherent) entangled state can be 
used for violation of the BI [19| . there are mixed (par- 
tially incoherent) entangled states which cannot violate 
the BI [20j . In fact, the ratio of the volume of the en- 
tangled states violating the BI to the volume of the en- 
tangled states obeying the BI is small. In particular, 
in the Hilbert-Schmidt metric this ratio equals 0.01085 
[l4j . This suggests that observation of the Bell inequal- 
ity violation is a more difficult task than observation of 
entanglement. Indeed, in the presence of decoherence 
the duration of the Bell-inequality violation is generally 
significantly shorter than the entanglement survival time 

B 

Until now most experiments on the BI violation have 
been performed with photons [2l], [22|, [23|, [2J] ■ However, 
recently there has been an increasing interest in testing 
the BI for various material systems, where decoherence is 
usually an important factor. In particular, experiments 
with ions in tra ps 251] , an atom-photon system [261 ] , and 
single neutrons [271 ] were performed. There are also var- 
ious theoretical proposals related to the BI violation in 
solid-state systems [H E3, IH l2l y0, [31j • Experiments 
on testing the BI in superconducting Josephson phase 
qubits [32| are currently underway [331 ]. 



Optimal experimental conditions for observation of the 
BI violation in superconducting phase qubits were con- 
sidered in Ref. 18. Both the ideal case and effects of 
various nonidealities, such as measurement errors and 
crosstalk [34], [35|], were analyzed in detail, while deco- 
herence was discussed briefly. In the recent paper [14[ 
the relation between entanglement and the BI violation 
was considered. In particular, the survival times of en- 
tanglement and the BI violation in the presence of lo- 
cal decoherence were studied, in the frame of a decoher- 
ence model more general than in previous publications 
H 0, S 0, I, S M EH, El El 111 E3, the focus be- 
ing on the class of "odd" two-qubit states relevant for 
experiments with superconducting phase qubits [34 , I36J ] . 
Note, however, that optimal detector configurations pro- 
viding the maximal BI violation were not discussed pre- 
viously, except for the case of pure dephasing [16[ (see 
also brief remarks in Refs. Ej and [181). The knowledge 
of all optimal configurations is important for planning 
experiments, since some detector configurations can be 
easier to realize than others [16j . 

In the present paper, we provide a comprehensive dis- 
cussion of effects of local decoherence on the BI violation. 
We obtain both optimal states and all detector configu- 
rations which yield the maximal BI violation. Our de- 
coherence model includes energy relaxation at the zero 
temperature (known also as spontaneous decay or am- 
plitude damping) and dephasing (phase damping). We 
analyze analytically and numerically the general case and 
a number of important special cases. In particular, we 
study the experimentally relevant classes of the "even" 
and "odd" states. We also discuss the combined effect of 
decoherence and local errors, basing on the treatment of 
errors in Ref. [l8l . 

The paper is organized as follows. In Sec. [TT]we discuss 
the BI and properties of maximally entangled states. In 
Sec. IIIII we obtain all optimal configurations of the de- 
tectors which maximize the BI violation for any given 
(pure or mixed) state. Section ITVl is devoted to effects of 
local (independent) decoherence of the qubits on the BI 
violation. In Sec. [V]we consider combined effects of er- 
rors and decoherence. Section lVll provides the concluding 
remarks. 



II. PRELIMINARIES 

In this section we introduce the notation and review 
properties of the BI and maximally entangled states. 
Here and below till Sec. [V] we assume the absence of 
measurement errors. 



A. The Bell inequality 

We consider a pair of qubits, i.e., two-level systems a 
and b. Each qubit has the states |0) and |1). Measure- 
ments of two qubits satisfy the Bell inequality, provided 
that a realistic (classical) theory holds and there is no 
communication between the qubits [15| . We consider the 
Clauser-Horne-Shimony-Holt version of the BI [37], l38j 



2< 5<2, 



(1) 



where 



S = E{a, b) - E(a, b') + E(a', b) + E(a' , b'). (2) 

Here a and a' (b and b') are the unit radius- vectors on the 
Bloch sphere along the measurement (detector) axes for 
qubit a (b) and E(a, b) is a correlator of the measurement 
results for the qubits a and b in configurations a and b, 
respectively, 

E(a,b) =pn(a,b) + p 00 (a, b) -p w (d 1 b) -p 01 (a,b), (3) 

where py (a, b) (i, j — 0, 1) is the joint probability of mea- 
suring results i and j for qubits a and 6, respectively. 

Note that the minus sign can be moved to any term in 
Eq. @ • This is related to the freedom of the interchange 
of the detectors for each qubit. Indeed, the substitution 
a «-> a' results in the interchange of the signs of the second 
and fourth terms in Eq. @. Similarly, the substitution 
b «-> b' (a «-> a', b <-> b') is equivalent to moving the minus 
sign to the first (third) term in Eq. ^. 



B. The Bell operator 

According to quantum mechanics in the absence of 
measurement errors 



PlJ (a,b) = T r {P?(a)PUb)p}. 



(4) 



Here p is the density matrix for the two qubits, P" = 
Pi® I 2l and P- — I 2 ® Pi, where /„ is the n x n identity 
matrix and Pi are the single-qubit projection operators, 

P 1 (a) = ^(I 2 + a-3), P (a) = ~(/ 2 - a ■ 3) , (5) 

3 = (cr x ,(Ty,<j z ) being the vector of the Pauli matrices 
[3|. Here and below we assume that |1) and |0) are the 
eigenvectors of a z with the eigenvalues 1 and — 1, respec- 
tively. 



Equations (O and JU) yield that E(a, b) = Tr(ABp). 
Here 

A = P«(a)-PS{a)=a-a a (6) 

and similarly B = b-3\,, where a a = 3®I 2 and a\, = I 2 (&3, 
the eigenvalues of A and B being ±1. Correspondingly, 
as follows from Eq. @, 



S = Tr(Bp), 
where the Bell operator [39] B is 

B = AB-AB' + A'B + A'B'. 



(7) 



(8) 



Here A 1 = a! -a a and B' = b' -3^. Note that measurements 
along a given direction a described by the projection op- 
erators ([5]) can be performed by means of a suitable ro- 
tation of the qubit and a subsequent measurement in the 
basis {|0),|1)} [HHHt 

In this paper we look for experimental conditions which 
are the most favorable for an observation of the BI viola- 
tion. Such conditions are reached when 15*1 — 2 is positive 
and maximal. We will use the following properties [18| 
of S which follow from Eqs. $7§ and (JU). 

(i) S is invariant under arbitrary local unitary trans- 
formations of the qubits, 



P -» (U a ® U b )p{Ul ® U, 



and the corresponding rotations of the detectors, 



(9a) 



R„a, 



R a a', b^R b b, b'^R b b', (9b) 



where U a (U b ) is a unitary matrix for qubit a (b) and R a 
(R b ) is the rotation matrix corresponding to U a {U b ), so 
that, e.g., U a (r a ■ 3)U\ = (R a r a ) ■ 3. A rotation matrix 
R is an orthogonal matrix, R T R — I3, with det(i?) = 1. 
(ii) S inverts the sign if the vectors a and a' (or b and 
b') invert the sign (e.g., a — > —a, a' — * —a!). Therefore, 
for a given state the maximal and minimal values of S 
are equal by the magnitude, yielding equal violations of 
the both bounds in the BI |T]). As a result, it is sufficient 
to discuss only the conditions for achieving the maximum 
of S. Below we denote by S+ and S- the maximum and 
minimum, respectively, of S for a given state and by 5 max 
the maximum of S over all states and detector axes. 



C. Maximally entangled states 

In the ideal case, when there is no decoherence or er- 
rors, the maximal and minimal values, S'max and Smm, 
which S can achieve are 14 lH 



2>/2, Smin = -2V2. 



(10) 



These limits are obtained for any maximally entangled 
state [43| . The BI violations are often considered for the 



following maximally entangled states, called also the Bell 
states 3, 



|* ± ) = (|11)±|00))/V2, 
|tt±) = (|10)±|01»/V2. 



(11) 
(12) 

For each maximally entangled state there are infinitely 
many optimal (i.e., producing a maximal BI violation) 
configurations of the detector axes a, a', &, and b'; all 
such configurations were described in Ref. UJJ In Sec. 
IIIII we describe the configurations maximizing S for an 
arbitrary (pure or mixed) state. 

It is useful to have a general expression for maximally 
entangled states in the "standard" basis of the qubit pair. 
It can be shown by performing the Schmidt decomposi- 
tion of the general two-qubit state that, with the ac- 
curacy to an overall phase, the most general form of a 
maximally entangled state is [43], E3 

|*) = d|ll) + c 2 e Ml |10) - c 2 e* Q2 |01) + cie l ( Ql+Q2 )|00) 

( 13 ) 
Here a\ and a 2 are real numbers, whereas ci, c 2 > and 

cf+cl = l/2. 

With the help of local rotations of the qubits around 

the z axis, Eq. (fT3"]) can be reduced to one of the two 

wavefunctions 



|*) =c*|$±) +c*|^ 



T/> 



(14) 



where either the upper or the lower signs should be used 
simultaneously, c$ and c* are any real numbers satisfying 
1, and |$±) and |<3>±) are the Bell states (JTrj 



and H2). Equation ([13]) is used in Sec. QV 



III. CONDITIONS FOR THE MAXIMAL BI 
VIOLATION IN A GIVEN STATE 

While formulas for the maximum S+ of S are known 
for both pure [45], El| and mixed [43] states, only one 
optimal detector configuration (i.e., a configuration for 
which S+ is realized) was provided in Refs. EH EflEZt In 
contrast, Samuelsson et al. [l|| showed that for the Bell 
state |$+) in the presence of dephasing there is a family 
of optimal configurations depending on one continuous 
parameter. In this section we extend the method of Ref. 
[43 in order to obtain all optimal detector configurations 
for any (pure or mixed) state. We show that the set 
of optimal detector configurations generally depends on 
one continuous and one discrete parameters, though in 
special cases the number of continuous parameters can 
equal two or three. 



A. Real representation of a two-qubit state 

It is useful to consider the following representation [47], 
EH, E2l of the two-qubit density matrix 

p = (h + r a ■ 5 a + n ■ 5 b + 5 a T5 b ) /4. (15) 



Here fk is the Bloch vector characterizing the reduced 
density matrix for the qubit k, so that, e.g., p a = Tr b p = 
{1 2 + r a ■ 5)/ 2, T is a matrix with the real elements 

T mn = Tr(per r a „crk ) (to, n = x,y, z), (16) 

and 5 a T5 b = J2 Z m ,n=x 7 ™™ <J m a ni where 5 k = 
(o*,o*j,o%) (k = a,b). 

Consider some useful properties of T. Since the eigen- 
values of the Hermitian operators cr^cr^ equal ±1, Eq. 
(TTB1 implies that 

\r mn \ < 1. (17) 

There exists the polar decomposition Q 



T = VVU, 



(18) 



where V is a 3 x 3 orthogonal matrix, V V = ^3, and 
U = T T T is a real symmetric matrix with nonnegative 
eigenvalues u\, u 2 , and 1*3; 113 being the smallest eigen- 
value (0 < U3 < ui, m 2 ). V is determined uniquely and 
given by V = TU~ 1 ' 2 , only if u 3 7^ 0; this is the most 
interesting case, as shown in Sec. IIIIB1 The determinant 
of V equals 1 or -1 for dct(T) > and det(T) < 0, re- 
spectively; when det(T) = (which means that u 3 = 0), 
V can be chosen such that det(V) = 1. 

Under a local unitary transformation p — > (U a ® 
Ub)p(Ul ® Ul) Eq. (TTSl) changes so that r k — ► R^k 
and 1481 



T 



R a TR b , 



(19) 



where R a ,b are defined after Eq. (|9bp. As follows from 
Eqs. (fT5)) and (flU)) . with suitable rotations of the qubits, 
R' a and R' b , the matrix T can be reduced to one of the 
two diagonal forms, T' = R' a T{R' b ) T = ±VU 7 - Here R' b is 
such that W = R' b U(R' b ) T is diagonal and R' a = ±R b V T . 
The plus and minus signs in the above formulas are ob- 
tained for det(T) > and det(T) < 0, respectively [on 
choosing det(V) = 1 when det(T) = 0]. As a conse- 
quence, in view of Eq. (Q~7]), we obtain 



< W3 < Ml, u 2 < 1. 



(20) 



Note that an arbitrary two-qubit state reduces to a 
simpler form by a local unitary transformation which di- 
agonalizes T. Inserting a general diagonal T into Eq. 
([j"5]) . we obtain that all states with a diagonal T have 
the form 



P11 P12 P13 Pu \ 

P21 P22 P23 P13 

P31 P23 P33 Pl2 

Pu P31 P21 Pa 



(21) 



where pij subscript values i,j = 1,2,3,4 correspond to 
the basis {|11), 1 10) , |01), |00)}. In the state ([2D_pi 2 = 
P34 and P13 = P24, whereas pu and p 2 3 are real [50(. 
In view of Eq. (fl6|) , for the state |2l|) we obtain 



T = 2 diag(p 23 + p u , P23 - Pu, 1/2 - P22 - ^33)- (22) 
Note that this expression is independent of pi 2 and P13. 



Maximal BI violation 



Let us review the derivation [47) of the maximum S+ 
of S for a given state. Inserting Eq. J5]| into ([7]) and 
taking into account Eqs. © and (|16p . we obtain that 



S = aTb- aTb' + a'Tb + a'Tb'. 



(23) 



The vectors 6 and b' can be always written in the form 



6 = 
6' 



c 



•; cos(C b /2) 4 
?i cos(C b /2) - 



?' 2 sin(C 6 /2), 
c' 2 sin(C 6 /2), 



(24) 



where c[ and c' 2 are orthogonal unit vectors and (& is the 
angle between b and 6' (0 < Cb < 7r). Inserting Eq. ([24]) 
into (H2]) yields 5 = 2[aTc 2 sin(Cb/2) + a'Tc[ cos(Cb/2)]. 
To maximize this expression, one should require a and a' 
to be parallel to Tc 2 and Tc[ , respectively [51| , yielding 



= Tc'/\Tc' 



i = Tc' 1 /\Tc' 1 \ 



Then maximizing S over £& results in 

C& = 2arctan(|Tc 2 |/|Tc' 1 |), 



(25) 



(26) 



where \v\ denotes the length of a vector v, and S = 

The maxi- 



■M^2- 



\Tc' 2 \ 2 = 2t/F\UF[ 
htained when (s 

c'iUci + c' 2 Uc' 2 = u\ + u 2 



2 v /|Tc' 1 | 2 + \-j Cgl' = -iy/c\uc\ +c 
mum of 5 is obtained when (see Sec. IIIIC|) 



yielding [47 



5 4 



2\A«i + u 2- 



Hence, the BI violation, S + > 2, occurs when ui + 
1. Equations ([2"0)) and ([2"8"| imply the limits on S+, 

0<S+< 2\/2. 



(27) 



(28) 



u 2 > 



(29) 



Here the lower limit, S+ = 0, is obtained for the states 
with T = 0, which are, in view of Eq. (|22| . the states 
locally equivalent to p — diag(/9ii,p22jl/2 — P22, 1/2 — 
/On). [In this case the elements pi 2 and P13 in Eq. |2T|) 
can be set to zero by local rotations of the vectors r a and 
r*h in Eq. (fT"5|) to the 2 axis.] Such states can be shown 
to be mixtures of product states with one of the qubits 
in the maximally mixed state I2/2. 

As an example, let us apply Eq. (|2"5)) to the state (f2"Tj) . 
In view of Eq. ((22j), U = T 2 , yielding 



2[max{8p^3 + 8pf 4) 4(|/t)23| 

+ (l-2/ )22 -2 /033 ) 2 }] 1/2 . 



|Pi 4 |) S 



(30) 



The states with det(T) > do not violate the BI. 
To show this, it is sufficient to consider a diagonal 
T = \/U, since, as mentioned in Sec. IIII A[ such T 
can be obtained for any state with det(T) > by 



means of local unitary transformations, which do not 
change S. A diagonal T = \fu is given by Eq. (|22|) 
with nonnegative matrix elements, which implies that 
r = P22 + P33 < 1/2 and P23 > \pii\- As follows from 
Eq. (dl|), TrU = TrT 2 = 4p 2 4 + 4p^ 3 + (1 - 2r) 2 . For 
a given r, TrU is maximal if p 2 3 and \pu\ are maximal 
under the above constraints, i.e., if \pu\ = P23 = r/2, 
which yields TrU < 6r 2 - 4r + 1 (0 < r < 1/2). This 
expression achieves the maximum TrU = 1 for r = 0. 
Thus, in the case det(T) > we have TrU < 1, the 
value TrU = 1 being obtained for the states which, with 
the accuracy to local unitary transformations, have the 
form p — diag(pn, 0,0,^44) and hence have T = U = 
diag(0,0, 1) [see Eq. $T2§}. In view of Eq. (J28J) , these 
states yield S+ — 2; they include, in particular, pure 
nonentangled states. 

However, the BI violation implies TrU > iti + U2 > 
1. Hence, a necessary condition for the BI violation is 
det(T) < 0. As a consequence, in view of Eq. (Tj"5)) . for 
the states violating the BI all u, do not vanish. 



C. Optimal detector configurations 

Consider detector configurations providing the maxi- 
mal S ([28]). Equation $2J§ obviously holds for c\ = ci 
and c 2 = c'2, where Ci (i — 1,2,3) are the unit orthog- 
onal eigenvectors of U corresponding to the eigenvalues 
Wj. This is the choice of Cj and c 2 which was made in 
Ref. [43 (similar choices were made also in Refs. [45] and 
1461 ). In this case Eq. (|2"6"1) yields Cfc = Co ; where 



Co = 2 arctan \Zu 2 /ui, 



(31) 



since \Tci\ 2 = cJJci — UiCi ■ c% = Uj. Hence, Eq. ([24 
becomes 



b = ci cos(Co/2) + c 2 sin(Co/2), 
b' = ci cos(C /2) - c 2 sin(Co/2), 

and, in view of Eqs. ([23)1 and ([X5]). 

a = e 2 , a' = e ! , 



where 



Tc,/ 



Vc, (* = 1,2). 



(32) 



(33) 



(34) 



The vectors ei^ are orthonormal, since V is an orthogonal 
matrix. 

The optimal detector configuration given by Eqs. ([3"2"[1 
and ([33]) is not unique. To obtain all possible de- 
tector configurations providing Eq. (|28p . we consider 
the derivation of the maximum (l27l). as follows. If ci 



is a unit vector orthogonal to c\ and c 2 , then c^ = 
J2j=iWijCj (i — 1,2,3), where W is an orthogonal 
3x3 matrix, W T W = I 3 . We have J^Li S W i = 



Sj=i CjMcj = Mi + ^2 + U3. Moreover, 



c 3 Wc 3 = W 3 >i + W 3 2 2 m 2 + W 3 2 3 u 3 > «3, (35) 

'- _i_v\)2 _i_v\;2 _ ^ and U3 < Ul + U2 _ As a result, 
c i^c j + c ' 2 Uc 2 = U1+U2+U3 — S' 3 US' 3 < mi + M2- Hence, 



since W% 1 + W 32 



the maximum (|27p is achieved when the expression 
is minimal, which occurs for W 3 i = W 3 2 = and W 33 = 
±1, i.e., for c 3 = ±c 3 . In this case c[ and c 2 are an 
arbitrary pair of orthonormal vectors in the plane defined 
by c\ and ci- All such c' x and c 2 are given by 



ci cos 77 ± c*2 sin 77, 



-ci sin 77 ± C2COS77. (36) 



Here 77 is the arbitrary angle of rotation of c\ and c 2 , and 
the two signs before c 2 correspond to a reflection with 
respect to the c\ axis. The above derivation implies that 
no pair of unit orthogonal vectors c\ and c 2 other than 
those in Eq. ([56]) can satisfy the condition ([2"T]) . unless 
m 3 equals u\ or U2 (the latter cases are discussed in Sec. 

IhTfI . 

Inserting Eq. ([36]) into Eqs. ([24]) and ([25]) and taking 
into account Eq. (|34[) . we obtain 

a= (— eiv^ sin 77 ± e-i^fui cos 77)/ 1 Tc 2 1, 

a' = (ei^tTcosr/i e 2 \/u2 sin 77)/ J Tcj |, 

6 = ci cos[?y + (b{v)/ 2 } ± c 2 sin[?7 + C&(»?)/2], 

b' = c x cos[?7 - C b (r])/2] ± c 2 sin[?7 - Cb(v)/% (37) 

In Eq. ([37]) the upper (or lower) signs should be used 
simultaneously. The quantities |'Z - c*^ | and \Tc' 2 \ are given 
by 



|Tc" 1(2) | = yj[ui +u 2 ± (mi - u 2 )cos2?7]/2, (38) 

where the upper (lower) sign corresponds to ITc'jJ 
(|Tc 2 |). Inserting Eq. ([38]) into Eq. ([26]) and performing 
trigonometric calculations yields 



&fa) 



arccos 



ui - u 2 

«1 + U-2 



cos2?7 



(39) 



Equation (|39|) can be compared with the angle £ a be- 
tween a and a', satisfying cos^ a — a- a'. As follows from 
Eqs. ([3"T]) and ([55]) and some calculations, 



Ca(v) = arccot 



U 2 - Ml 



2^/uTu2 



sin 2r; 



(40) 



(we assume that < Ca,0> < 7r )- Equations ([39]) and 
([40]) imply that the angles ( a and £& vary with 77 between 
the values Co (J3T]) and 7r — Co with the period 7r. In 
particular, for 77 = 0,±7r/2,±7r (77 = ±7r/4, ±37r/4) one 
gets ( a — tt/2 (C& = tt/2), whereas (b {(a) acquires a 
maximal or minimal value. Moreover, Eqs. (|37[) - ([3T)]) 
imply the relations 

a?(r) ± tt/2) = ±a( V ), b'(r) ± tt/2) = ±b(r)), (41) 



whereas a, a', b, and 6' change the sign for 77 — > 77 + n. 

Note that the optimal detector orientations ([3"7| de- 
pend on mi and it 2 only through the ratio U2/U1. As 
a result, two different states with the same V, ci, C2, 
and 112/ui have the same or, at least, overlapping sets of 
optimal detector configurations. (The overlap may be in- 
complete only when, at least, for one of the states u\ = u 3 
or u 2 = u 3 , see Sec. IIIIFI ) In particular, for the states 
p and p 1 with the matrices T and T', respectively, satis- 
fying T' = fT (/ > 0), the respective quantities S + and 
S' + obey S' + = fS+, whereas the optimal configurations 
are the same for the both states. An example of such 
states p and p' are, respectively, the input and output 
states of the depolarizing channel Q , 



p' = fp+(l-f)(h/4) (0 </<!). 



D. Polar coordinates 



(42) 



According to Eq. ([3~T]) , the optimal detector directions 
for qubits a and b are confined to the planes (e\, e^) and 
(ci,C2), respectively. It is convenient to specify these 
directions by means of polar angles. 

To this end, we introduce the polar coordinate v in the 
(ci,C2) plane, which is counted from c\ in the direction 
where v = 7r/2 corresponds to C2, and the polar coordi- 
nate (5 in the (ei,e*2) plane [differing generally from the 
(ci,C2) plane], which is counted from e\ in the direction 
where S — tt/2 is the polar coordinate of €2- Then Eq. 
([37]) can be recast in the form 

a = e*i cos[S a (r])} ± e 2 sin[<5 a (77)], 

a! = ei cosK(?7)] ± e 2 sin^ry)], 

b — c 1 cos[p b (r])} ± c 2 sm[v b (i])], 

V = Si cosK(7?)] ± c 2 sm[iy' b {r,)}. (43) 

Here S a , 5' a , v b , and v' h {-8 a , —S' a , -v b , and —v' b ) are 
the polar coordinates of the vectors a, a', b, and b', re- 
spectively, when the upper (lower) sign in Eq. (|43[) is 
realized. The configurations corresponding to the two 
different choices of the sign in Eq. ([37]) or (|43|) trans- 
form to each other by the reflection of the detector axes 
for qubits a and b with respect to the axes e\ and ci, 
respectively. 

As follows from Eq. (J4"T]) . the functions 

S a (rj), S' a (rj), 7^(77), and ^(77) satisfy the relations 

6 a (v) = S' a (v + rr/2), ^(77) = i/ b (r, + tt/2). (44) 

Expressions for these functions are obtained by compar- 
ing Eqs. ([37]) and ([43]) . In particular, we obtain 



(45) 



Mv) = V + CbW/2, v' b (v) =V~ Ur])/2. 



This implies that 77 = (vb + l/ ' b )/2, yielding a geometric 
interpretation of 77: 77 equals the polar coordinate of the 
bisector of the angle between b and b' when the upper 
sign in Eq. ([57]) or ([4"3"]l is realized. 



Let us now obtain S' a (r]); then 8 a (rf) will follow by 
the first Eq. [[31]), The second Eq. (J37J) implies that 
cos S' a = a' • e*i = y^cos^/ITc^l, whereas the sign of 
sin<5^ coincides with the sign of sin??. With the account 
of Eq. (|38]), this yields (for -tt/2 < 77 < vr/2) 

^(77) = arctan f y / u 2 /u 1 tan 77 J , ^(77) = 5' a (rj + ir/2), 

(46) 
It is convenient to consider the polar angles as continuous 
functions of 77 without restricting them to any interval. 
In particular, S' a (ri) in Eq. (|46p can be extended continu- 
ously beyond the interval — tt/2 < 77 < tt/2 with the help 
of the equality S' a (r] + kn) = 5' a (n) + kTT (k = ±1,±2, ...)• 
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FIG. 1: Polar angles of the detector axes versus r\ for 112/ui = 
0.3 (Co ~ 1-00), as given by Eqs. fl45j) and fl46j). 

Figure [1] shows the dependence on 77 of the polar co- 
ordinates of the observation axes when the upper sign is 
chosen in Eq. J43J and uiju\ = 0.3 (Co ~ 1-00). Note 
that in Fig. Q]the observation angles are not restricted to 
an interval of the length 2tt to avoid discontinuities in the 
plots. Our calculations show that the set of the polar co- 
ordinates of the detectors taken as continuous functions 
of v is always ordered as follows (cf. Fig. Q]), 



8 a > Vb > S' a > v' b . 



(47) 



E. 



Symmetry with respect to the exchange of the 
qubits 



In the above derivation (Sec. IIII F3|) the detectors for 
qubit b are treated differently than those for qubit a [cf. 
Eqs. ([2"4"]l and (|25p], resulting in apparently different so- 
lutions for the two qubits [cf. Eqs. ([31]) and (|3tj])]. It is of 
interest to check whether these solutions are symmetric 
with respect to the qubit swap. We cast Eq. (|23[) as 



As follows from Eq. (18]), brT a = bV T TV T a 
{Vb)T(V T a). Thus, Eq. (gSJ) can be recast as 



S = (Vb)T(V T a) - (Vb')T(V T a) 
+{Vb)T{V T a') + (Vb')T(V T a') 



(49) 



Comparing Eqs. ([23]) and ([49]) , we obtain that S is in- 
variant under the simultaneous substitutions 



Vb' 



a, 



Vb 



V T a' 



V T a 



b'. 



(50) 



In view of Eq. (|50p . we obtain that the qubit-swap 
operation is equivalent to the replacement 

8 a -> v'b, S'a -> v b , v h -v 8' a , v' b -* 8 a (51) 



in Eq. (gSJ). As shown in Sec. |TTCl Eq. ([57 jl or (|^5 ]) de- 
scribes all optimal detector configurations (at least, for 
""3 < U\,U2)- Therefore, optimal configurations trans- 
formed by the replacement (|5ip should belong to the 
set of detector configurations described by Eq. (f3"T)) or 
(I43p . This is confirmed by our numerical calculations 
which show that 8 a , 8' a , and v b as functions of v' h coin- 
cide, respectively, with — p' b , —v , and — 8' a as functions 
of —8 a . In other words, Eqs. (|4"5|) and (|46[) imply that 
if (6 a ,8' a ,Vb,v' b ) are optimal polar angles of the vectors 
a, a', b, and b' , respectively, then (— u b , — Z/&, — 8' a , — 8 a ) 
[and hence (i/ b , v , S' a , S a )] are also optimal polar angles 
for the respective vectors. Thus, the set of the detector 
configurations (|57|) or (|4T{|) is invariant with respect to 
the qubit-swap symmetry relation (|50p [or (|5ipj. 

As follows from the above, the set of all optimal config- 
urations maximizing S for a given state generally depends 
on one continuous parameter (77) and one discrete param- 
eter [which corresponds to the two possible signs in Eq. 
(|3T|) or ([32])] ■ However, when 7^3 is equal to u\ and/or u 2 
(a degenerate case), the set of optimal configurations is 
characterized by two or three continuous parameters, as 
discussed in Sec IIII Fl 



F. Degenerate cases 

Consider the optimal detector configurations in degen- 
erate cases. 

a. The singlet state. We begin with the singlet state 
p = |W_)(^_| [Eq. p2p ]. which is often considered in 
connection with the BI. In this case Eq. (fll)]) yields an 
especially simple result, T = —I 3 and U = T 2 = 7 3 , 
implying u\ = u 2 = U3 = 1. Hence, in view of Eq. (|28p . 
S+ = 2\/2 [cf. Eq. pop]. Now any two perpendicular 
unit vectors can be chosen as c\ and c 2 . Since V = T = 
-I 3 [see Eq. ([15])] and Co = tt/2 [see Eq. ([3T])]. in view 
of Eqs. ([55]) and ([3^]) . we obtain a = — c 2 , a' = — C\, and 



b = -(a + a')/V2, b' = (a-a')/V2. 



(52) 



bT 1 



bT 1 



b'T 1 



b'T T c 



(48) 



Equation (|52p with arbitrary mutually perpendicular a 
and a' is known to provide all optimal configurations for 



the singlet state [l8|. They are characterized by three 
continuous parameters, since all optimal configurations 
can be obtained by arbitrary rotations of any given opti- 
mal configuration [18| . 

b. Maximally entangled states. Any maximally en- 
tangled state is obtained from the singlet by local ro- 
tations of the qubits, so that, in view of Eq. (fT§|) , for 
such a state T = — R, where R = R a Rj corresponds 
to the rotation of qubit a relative to qubit b. In this 
case r a — f% = 0, and hence, in view of Eq. (|15[) . 
p = (J 4 — <j a R<Jtj)/A, i.e., a maximally entangled state 
is determined [44[ by R. 

Now U = R T R = I3, as for the singlet. Proceeding 
as in the case of the singlet, we obtain a = —Rc%, a' = 
—Rc\. As a result, all optimal configurations providing 
S + = 2\/2 are given by Eq. (f5"2")l with a — > R T a and 
a! — > R T a' , so that 

b= -(R T a + R T a')/V2, b' = (R T a-R T a')/V2. (53) 

Since the rotation matrix i? is determined by three pa- 
rameters (e.g., the Euler angles), all optimal configura- 
tions for all maximally entangled states are characterized 
by six continuous parameters [18J ]. 

Note in passing that all configurations producing S- = 
— 2\/2 follow from Eq. ((53|) with a — ► —a and a' — ► —a' 
(see Sec. Ill B|) . yielding 

b = {R T a + R T a')/V2, b' = (R T a' - R T a)/V2. (54) 

c. Completely degenerate case (u± = U2 = W3J. Let 
W = u/3, yielding ui = U2 = U3 = u. In this case Eq. 

yields S + = 2>/2u, so that the BI is violated for 



u > 1/2. 



(55) 



Now T = y/uV, where V is an orthogonal matrix [see 
Eq. fig])] . When the BI is violated [and hence Eq. l[55j) 
holds], det(T) < and hence i? = — V is a rotation ma- 
trix. In this case the optimal configurations are the same 
as for the maximally entangled state, Eq. ([53")) (cf. the 
last paragraph of Sec. IIIIC|) . As a result, in the com- 
pletely degenerate case the set of optimal configurations 
is characterized by three continuous parameters. 

An example of a state corresponding to the completely 
degenerate case is the Werner state [2(| 



p = V«|tt-><*-l + (l-V«)(l4/4), 



(56) 



which is a special case of the state (|4"2"j) . The BI violation 
condition (|55"|) for this state was obtained in Ref. |47I . 
For the state (|5rj|) T = — \ful$, yielding R = I3, and we 
obtain the same optimal configurations (|52|) as for the 
singlet. 

cf. Case u\ = U3 or «2 = U3. In this case we as- 
sume for definiteness that u-i =113. Then c\ is defined 
uniquely, whereas 02 can be any unit radius- vector in the 
plane perpendicular to c\ . Now the set of optimal config- 
urations ([3T)l or (|4^|) is characterized by two continuous 
parameters, rj and an angle specifying the direction of Ci 



in the plane (02,03) with respect to some reference axis. 
[In this case the discrete parameter is superseded by the 
new continuous parameter; indeed, now one of the two 
possible signs in Eqs. (|3T| and (|4*3")> can be omitted, since 
it is recovered for ci — » — ci\ 

e. Case u\ = ui- In the case u\ = U2 = u Eq. (|28j) 
yields 5+ = 2>/2u, as in the completely degenerate case. 
Equations pS jl -pP ]) yield \Tc[\ = |Tc' 2 | = 01 and 



Ca = Cb = 7T/2. 



(57) 



Now any two perpendicular unit vectors can be chosen 
as c\ and c*2 in the plane spanned by the eigenvectors of 
hi corresponding to u\ and 1*2- Hence, as follows from 
Eqs. (f3"2"| and (|33[) . the optimal configurations are given 
by Eq. (|53[) . where i? = — V and a and a' are arbitrary 
perpendicular unit vectors in the plane (e\, e-i). In other 
words, now in the case of interest det(V) = —1 the opti- 
mal configurations are a subset of the set of the optimal 
configurations for the maximally entangled state char- 
acterized by the rotation matrix R = — V. The optimal 
configurations are described by the following polar angles 
[see Eqs. (g5j) and 



5 a = V + tt/2, 5' a = 1-1, iy b = ij + tt/4, i/ h =n~ tt/4. (58) 

Now the set of optimal configurations is characterized by 
one continuous parameter and one discrete parameter as 
in the nondegenerate case. 



G. Special cases 

Here we consider the states with T assuming one of 
the two simple forms, 



T = <&&g(T x ,T x ,-T z ), 
T = dl&g(T x ,-T x ,T z ). 



(59) 
(60) 



In view of Eq. (j22|) . the most general two-qubit states 
with T given by Eq. (|59")) or (|6"0|) are described by Eq. 
(l2Tj) with P14 = or P23 = 0, respectively, so that in 
Eq. dM]) [dliD] t x = 2p 23 and r z = 2p 22 + 2p 33 - 1 
(tjc = 2pi4 and r 2 = 1 — 2p22 — 2/933). We assume below 
the validity of a necessary condition for the BI violation, 
det(T) < (Sec. IIIIBI) . which implies t x ^ and t z > 
0. The latter inequality is equivalent to the conditions 
P22 + P33 > 1/2 and p 2 2 + P33 < 1/2 for the cases (JS5J) 
and (|60|) . respectively. Without a loss of generality, we 
focus on the states with t x > in Eqs. ([55]) and (|6"0|) 
0. 

The cases ([591 and (pDl are of interest by themselves 
and are also of relevance below. In the both cases U = 
diag(rj , tJ , r|) , yielding, in view of Eq. 



(61) 



5+ = 2 max U/2t x , \/ t x+ t z} , 



i.e., S + = 2\[2t x for r^ > r z and 5+ = 2-^/tJ + t| for 
Tx < T"z- However, the optimal detector configurations 
are different in the cases (J5TJ)) and (|6^|) . as follows. 



Case (M 



In the case 



V = diag(l,l,-l). If t x > t z 



u\ = u 2 



tz., and we can choose C\ = x and c 2 



y, 



where x, y, and z denote the unit vector along the corre- 
sponding axis. Now ei (e^) coincides with c\ (62), yield- 
ing 5 = v = (j), where <fi is the polar coordinate in the 
horizontal (xy) plane. The optimal configurations lie in 
the horizontal plane, being given by 



(62) 



, <j>' a , <j> b , ft) = ±(0, tt/2, tt/4, 3tt/4) + C, 



where C is an arbitrary real number. Equation (|62p fol- 
lows from Eq. (|58|) . where we introduced the double sign 
to take into account the both signs in Eq. ([43]) . Hence, 
Eq. (|62|) describes all optimal configurations in the hor- 
izontal plane. 

If t x <t z , we choose 



Ml 



u 2 = u 3 = t x . 



(63) 



Then c*i = — e\ = z. Since u 2 — u 3 (see Sec. IIIIFI d). 
we choose c 2 = e 2 — xcos^o + ysintfio, where <f>Q is an 
arbitrary number which equals the polar angle of c 2 in 
the xy plane. The above expressions for c*x, c 2 , e\, and 
e 2 imply 

v = t: - 5 = 6, (64) 

where 9 is the polar angle in a vertical plane (a plane 
passing through the z axis); is counted from the z axis 
in the direction where 8 = tt/2 corresponds to c 2 . The 
polar coordinates of the optimal detector axes 8 a , 6' a , 9b, 
and Q' b as functions of the parameter 77 can be obtained 
from Eqs. @5J) 5 flUD, flB3J}, and ((Ml)- Thus, now op- 
timal configurations lie in any vertical plane; they are 
characterized by two continuous parameters, r\ and 4>o, 
in agreement with Sec. IIIIFI d. 

As discussed above, the optimal configurations are hor- 
izontal for t x > t z and vertical for t x < t z . For t x = t z 
we have the completely degenerate case (Sec. IIIIFI c) 
with R = — V = diag(— 1, — 1, 1) which corresponds to 
the rotation by n around the z axis. This rotation of 
qubit a yields \^-) —> \^+)- Hence, in this case the op- 
timal configurations are given by those for the Bell state 
|*+) <HU> (see Sec. glTFl b). In particular, Eq. §2$ pro- 
vides the horizontal optimal configurations [18( for |\&+), 
whereas the vertical optimal configurations for |vj/+) fol- 
low from Eqs. ([58]) and ([64]) and can be cast as [18| 



(6 a ,6' a ) = ±(0, tt/2) - C, (6 b , 9' b ) = ±(37r/4, tt/4) + C. 

(65) 
Here we introduced the double sign as in Eq. (|62"1) . 



2. Case [6$ 



In the case ([50]) V = diag(l, -1, 1). If t x > t z 



Ul 



~~ c*2 = y, yielding v = —5 = <f>. In this case the optimal 
configurations lie in the horizontal plane and are given 
by Eq. (B5]) with <\> a -> -<f> a and <j>' a -> -^, i.e., by 

(0 O> ^) = ±(0, -tt/2) - C, (^ 6 , ^) = ±(n/4, 3tt/4) +C. 

(66) 
For Ta, > r^ Eq. ([6"3"]> holds. As above, ci = ei = z and 
c*2 = x cos 0o + ysin ^o , but e*2 = x cos </>o — ysm <j>$ . Now 
the optimal detector axes for qubits a and b lie generally 
in different vertical planes characterized, respectively, by 
the polar coordinates <p a and <f>b (in the xy plane) of c 2 
and e*2, respectively, such that 

4>b = -0a = 00- (67) 

For a given value of 4>q, the optimal configurations are 
given by Eq. (|58[) . They are planar in two cases. For 
(J)q = and 7r the optimal configurations lie in the xz 
plane and satisfy 

v = 5 = 6, (68) 

whereas for <fio = tt/2 and 37r/2 they lie in the yz plane 
and satisfy 



v = -S = 6. 



(69) 



In Eqs. ([6"B"]) and (|69[) 6* is counted from z in the direction 
where = 7r/2 corresponds to x and y, respectively. 

In the intermediate case t x = t 2 the optimal configu- 
rations coincide with those for the maximally entangled 
state characterized by R = —V = diag(— 1, 1, — 1) (see 
Sec. IIIIFI bV This state, obtained from the singlet by 
the 7r rotation of qubit a around the y axis, is | < I > +) [Eq. 
(|11[) ] . In particular, the horizontal optimal configurations 
for |$+) are given by Eq. ([6"6"]) . The vertical optimal con- 
figurations are given by Eq. (I5S 



(70) 



(5 a ,6' a , v h , v' b ) = ±(0, tt/2, tt/4, 3tt/4) + C. 



"2 



t x , and we can choose c\ = e\ = x and c 2 



[The double sign is introduced here, as in Eq. ([rJ^]) . to 
describe all optimal configurations.] In the planar cases 
Eq. ([70]) simplifies according to Eqs. ([6"8"]l and ([6"0]) . 

Note that Eq. (JSOJ) follows from Eq. (HPJ under the 
7r rotation of qubit a around the x axis. As a result, 
the above optimal configurations for the case ([6"0]) can be 
obtained from those discussed for the case ([59"]) by the 
■k rotation of the detector axes for qubit a around the x 



3. Nonmaximally entangled pure states 

As a simple application of the above results, consider 
nonmaximally entangled pure states. Let us begin with 
the "odd" state 

|*) =cos/5|10)+sin/3|01) (0 < (3 < tt/2). (71) 

This state [as well as the states ([93]) and (|105[) below] 
is relevant for experiments with superconducting phase 



qubits \M,\3_, 
t x = sin 2/3 and t, 



For Eq. (n} T is given by Eq. (159 
- ' Now t x <t z , yielding [45j,|46] 



with 



= 1. 



S+ = 2A/l + sin^2/3. 



(72) 



For (3 = 7r/4 the state (|71"|) is maximally entangled and 
coincides with the Bell state \^>+) (see Sec. IIII G 1|) . For 
the nonmaximally entangled state ([7T|) (/3 ^ 7r/4) the op- 
timal configurations lie in a vertical plane and depend on 
two continuous parameters, as discussed in Sec. IIII G 11 
where now u\ = 1 and u 2 = sin 2/3. 

Any entangled pure two-qubit state can be obtained 
from the state (fTTj) by local rotations R a and Rb of 
the qubits a and &, respectively. These rotations satisfy 
V = RaTRj, where T and T = diag(sin2/3, sin 2/3, -1) 
are the matrices (fi~6| for the given state and Eq. (fTTj) . 
respectively. In view of the invariance of S under local 
rotations ([9]), we obtain that for an arbitrary entangled 
pure state, S + is given by Eq. ([72"]) , whereas the optimal 
configurations coincide with those for the state ([7T|) with 
the substitutions (|9b|) . 

In particular, consider the "even" state 

|$) =cos/3|ll)+sin/3|00) (0 < < tt/2). (73) 



It satisfies Eq. (|60l) and has the same r x , r z , and 5+ as 
the state (17T1) . For /3 = 7r/4 the state (|73|) coincides with 
|$+) (see Sec. IIIIG2JI , For /3 ^ tt/4 the optimal configu- 
rations are vertical and are obtained as discussed in Sec. 
IIII G 21 using u\ = 1 and u 2 = sin 2 2/3. These configu- 
rations can be obtained from the optimal configurations 
for the state ([7T]) by the 7r rotation of the detector axes 
for qubit a around the x axis (cf. the last paragraph in 
Sec. IIII G 21) . 



IV. EFFECTS OF DECOHERENCE 

A. Description of decoherence 

To investigate effects of decoherence on the BI viola- 
tion, we assume the following simplified picture of the 
experiment: after a fast preparation of the initial state 
Po, the qubits undergo decoherence during time t result- 
ing in the state p, then a fast measurement follows. Now 
in Eq. 



p = C(po) 



(74) 



where the superoperator (linear map) £ describes deco- 
herence of the qubit pair. We assume independent (local) 
decoherence of each qubit and the absence of any other 
evolution, so that 



transitions |1) — > |0)) and pure dephasing. The assump- 
tion of the zero temperature, T — 0, simplifies the for- 
mulas below. It is applicable to low-temperature systems 
(such as superconducting phase qubits), with ksT <C E q , 
where ks is the Boltzmann constant and E q is the en- 
ergy separation of the qubit. Moreover, the BI violation 
conditions were found in Ref. [lj to depend very weakly 
on the temperature. Under the above assumptions, the 
elements of the density matrix pk(t) of qubit k (k — a, b) 
obey the Bloch equations [52[ 



Pii(t) - -Poo(t) = -ptiW/Tf, 
Pi (t)=-Pio(t)/T 2 k , P k m (i) = 



-p k 01 (t)/T 2 k - (76) 



Here Tf and T 2 fc 

rpk 

L 2 



are the decoherence times, obeying 
< 2Ti , where the inequality occurs in the pres- 
ence of (pure) dephasing which proceeds with the rate 
r d = 1 / T 2 - 1 /( 2T i)- In the derivation of Eqs. J7B]) it is 
usually assumed that E q /Ti > 1/Tf, 1/T|. This condi- 
tion is satisfied in many systems. In particular, it holds 
for superconducting phase qubits. 

Equations (|76|) can be easily solved, providing the lin- 
ear map Pk{t) = £fe(pfe(0)). The superoperator £& can 
be written in terms of the Kraus operators (the operator- 
product form Q), 



C k (p k (0)) = J2K?p k {0)(KN, 



(77) 



where the Kraus operators are the special case for T = 
of the Kraus operators obtained in Ref. UA . 



Here 



Ik 



K\ 



Kl = 



-t/T£ 









s/T=7k~ 



o yi - ^ 



, Pk = Xk/y/lit 



K 



Tk o 

Pk 



-r5t 



Xk 



-tin 



(78) 



(79) 



The Kraus operators ([75)1 take into account energy re- 
laxation and dephasing simultaneously, extending thus 
the previously known Kraus operators |3fl which describe 
either effect separately. 

Combining Eqs. dJl, 1(75]). and d77J yields the expres- 
sion for the two-qubit superoperator 



£(po 



3 



KijPoK. 



Kij = K? 



>K 



(80) 



£ = £ a <g> C b 



(75) 



We consider Markovian decoherence which involves en- 
ergy relaxation at the zero temperature (i.e., spontaneous 



through the nine two-qubit Kraus operators Kij. As a 
result of decoherence, the initial two-qubit density matrix 
Po = {p'ij} evolves after time t to 



10 



lalbPn 7ahPl2 ^albPu 

la\P2l JailbPn +P22) X a X b p% 3 

AaTfoPsi X a X b p% 2 7b(7>?i+P3 3 ) 



A a A{,p 14 



7aP22 



7^3 



(81) 



p\, 



r 



where 7^/ b \ = 1 — 7a(b) ■ The ordering of the elements of 
the density matrix is specified after Eq. (J^TJ) . 

Decoherence generally breaks the invariance of S with 
respect to local transformations of po and the correspond- 
ing rotations of the detectors [cf. Eq. ©]. Aa a result, 
locally equivalent initial states may yield different maxi- 
mal violations of the BI. However, in the present model 
of decoherence S is still invariant under local rotations of 
the initial state around the z axis and the corresponding 
rotations of the detectors. 



However, due to the invariance of S with respect to ro- 
tations around the z axis, we can reduce the number of 
the fitting parameters from eight to six by setting, say, 
a y = by = 0. Combining the above invariance with the 
freedom of rotations of the detectors (see Sec. IIII1 es- 
pecially Fig. [1]) yields that one of the detectors can be 
fixed in any position, say, along the x axis, so that, e.g., 
a = (1, 0, 0). This further reduces the number of the fit- 
ting parameters to five. The second method involves the 
analytical approach of Sec. IIII1 as discussed in Sec. IIVCI 



B. Bell operator modified by decoherence 



C. States maximizing the BI violation 



It is useful to recast S = Tr (Bp) in the form S 
Tr(BdPo) with the modified Bell operator Bd = {C a 
C* b ){B) or 



B d = A d B d - A d B' d + A' d B d 



A' d B' d . 



(82) 



Here C k (fc = a, b) is the map adjoint (dual) to Ck that 
moves observables of the quantum system [53j and 



A„ 



C-l(A) 



{Kf^AKf, 











1= 


= 1 




etc. 


From Eqs. 


®, 


flU, 


and 


(|83| we 


obtain 


-4, ; 


= (7a-l)a z +9a-0 ; a 


) Qa 


= (A a Oa; 


, A a dy 7 



(83) 



where Xk is defined in Eq. (|79p . Expressions similar to 
Eq. (JH3I) hold also for A' d , B d , and B' d . The maximal vi- 
olation of the BI occurs always for a pure initial state (an 
eigenvector of the Hermitian operator B d corresponding 
to the maximal eigenvalue). 

We performed numerical calculations of the maximum 
S'max of S over all the states and observation directions as 
a function of the decoherence parameters with the help 
of the Mathematica 6 routine NMaximize. We used two 
methods. The first method is based on the fact that .Smax 
is equal to the maximum of the greatest eigenvalue of B d 
over the directions a, a', &, and b' , the optimal state be- 
ing given by the corresponding eigenvector. The detector 
axes are determined by eight independent parameters. 



Our numerical calculations by the method of Sec. IIVBI 
show that S'max can be always obtained for a pure state 
pa = |*)(*| of the form 



|*) = ci|ll) + ca|10) + C3|01) + C4|00), 



(85) 



where Ci are real coefficients such that J2i=i c f = 1- Note 
that the state l|55j) is not the most general two-qubit 
state, even if one takes into account the invariance of 
S with regard to local z rotations. Indeed, the general 
two-qubit state depends on six independent parameters, 
and local z rotations can decrease this number to four, 
whereas Eq. ([85]) depends on three parameters. 

If the state (|85p with certain values of c% provides SVax 
for some choice of the decoherence parameters, so does 
the state |*(a a ,a b )) = U z (a a ) <g> U z (a b )\^), where a a 
and otb are any real numbers and U z (a) = e ~ la,J z/ 2 ro- 
tates a qubit around the z axis by angle a. With the 
accuracy to an overall phase factor, 



\*(aa,a b ))= Cl |ll)+c 2 e Mb |10)+c 3 e^|01) 



J(a a +a b ) 



c 4 |00). 



(86) 



All optimal detector configurations for this state are ob- 
tained from those for the state ([55]) by rotating the de- 
tectors for qubits a and b around the z axis by the angles 
a a and a&, respectively. 

It is convenient for numerical calculations to express 
the coefficients c, through the parameters «i, K2, and K3 
by 



J 



(ci, C2, C3, C4) = (sinKi, coski sinK2,cosK;i COSK2 SUIK3, cos K\ cos k 2 cos K3). 



(87) 
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Now in Eq. (|8Tj) the elements p° = CiCj are real and we obtain from Eq. (p~6|) the matrix 



T 



KMpIz+Pu) 






AaA b (p°3 - Pu) 



\ a [p° 2i + (l-2 lb )p< 







A b [p° 4 + (l-2 7Q )p? 2 ] l/2-dA-7«/& 



76P33 



(88) 



r 



where d = 7 a + 7b ^ 27 a 7b (note that < d < 1). In the 
matrix T (|83|) only the xz and zx off-diagonal elements 
are nonvanishing. This allows us to obtain an analytical 
solution for S + and the optimal configurations, using the 
formalism of Sec. IIII1 as follows. 

From Eq. (|88[) we obtain that the nonzero elements of 
U = T T T equal 



U xx = 4A 2 {A 2 ( Pl ' 4 + p° 23 ) 2 + [p° 34 + (1 - 2 7q V? 2 ] 2 }, 

W w = 4A 2 A 2 (p i , 4 -p 2 ! 3 ) 2 , 

U zz = 4A 2 [p 2 l 4 + (l-2 7 b)p? 3 ] 2 + 5 2 , 

«** = Z4* = 4A 2 Ab(p? 4 + p° 3 )[p° 4 + (1 - 2 lb )p° 13 ] 

-2\ b [p° 34 + (l-2 7a ) p° 12 }g, (89) 

where g = 1 — 2dp\ 1 — 27 Q p22 — 27bP 33 . It is easy to see 
that the eigenvalues of U equal hi yy and 



U± = {U xx + U zz )/2 ± ^{U xx -U zz y/A + W xzl (90) 
the corresponding eigenvectors being y and 

c± = [(U± -U xx f +U 2 XZ ]- 1 ' 2 {U XZ ,Q,U ± -U XX ). (91) 
Thus, we obtain that 



til = Ua 



u 2 



m&x{U-,Uyy}, 



(92) 



whereas c\ = c+ and c 2 — C*_ if u 2 = U- or c 2 = y if 
u 2 = Uyy. These values of u\, u 2 , cf, and c 2 provide 
Sa- and the optimal detector configurations for the state 
(pj). on employing Eqs. Ip ]l . ip| ) . (pj) . and ((37])- (J39J) 
[or (O, (gSJ), and |15))]. Generally, the optimal detector 
axes lie in the xz plane, when c 2 = c_, or in two planes 
passing through the y axes, when c 2 — y (since then 

e*2 = y)- 

We use this solution to obtain S nmx and the corre- 
sponding optimal state, by maximizing 5+ numerically 
over the parameters «i, k 2 , and k 3 . This procedure is 
significantly faster than the numerical method described 
in Sec. IIVB1 Before the discussion of the results of nu- 
merical calculations in Sec. IIVE1 we consider important 
cases which admit simple analytical solutions. 



D. Analytical solutions 

1. Horizontal optimal configurations 

An especially simple solution is obtained when the op- 
timal observation axes lie in the xy (horizontal) plane. 



Then Eqs. ^ and (JUJ) with a z = yield B d = A a A b £, 
and hence S = X a XbSo, where So = Tr(Bpo) is obtained 
in the absence of decoherence (Sec. [II]). As a result, 
the value of S maximized over all states and horizontal 
observation directions is 



Sh = 2v 2X a Xf, 



(93) 



This value is obtained only for the maximally entan- 
gled states which have horizontal optimal detector con- 
figurations in the ideal case, as, e.g., the states |vl / +) and 
|$ + ) discussed in Sec. IIII Gl All such states are given by 
the expressions 



|*) = (|10)+ e M |01))/\/2, 
|$) = (|ll)+e M |00))/V2. 



(94) 
(95) 



Indeed, taking into account that S is invariant under 
identical rotations of qubits and detectors, Eq. ^j, and 
that all maximally entangled states are related by a ro- 
tation of one of the qubits (cf. Sec. IIII Fi b), we obtain 
that all maximally entangled states with horizontal opti- 
mal configurations result from the state |^+) (or 1$+)) 
on applying to one of the qubits an arbitrary rotation 
around the z axis and possibly a 7r rotation around the 
x axis, since only such rotations do not take the detector 
axes out of the horizontal plane. All the resulting states 
are given by Eqs. (j94|) and ([95]) . 

The states (|9"4"|) and (|9"5"|) are special cases of the states 
(I105P and (|106p . which are discussed in detail in Sec. 
IIVD2I In particular, in Sec. IIVD 21 the validity condi- 
tions of Eq. (|93| are obtained. 



2. Even and odd states 

It is of interest to consider maximal violations of the 
BI for the special cases of Eq. (1851) , the classes of "odd" 
and "even" nonmaximally entangled states, Eqs. ((7i"|) 
and (|T3")) . respectively. For the odd and even states, re- 
spectively, Eq. (|88p yields 



T = diag(A a Af, sin 2/3, A a Ab sin 2/3, 1 — 7+ 



7_ cos 2/3), 
(96) 



T = diag(A a A(,sin2/3, — A a AbSin2/3, 1 — d — d cos 2/3), 

(97) 
where j± = j a ± 7b- Equations (f9"o) and (|9"7| have 
the form of Eqs. (fS"9|) and (f60|) . respectively, with 
t x = A a Ab sin 2/3, whereas r 2 = 7+ + 7_ cos 2/3 — 1 and 
t z = 1 — d — dcos 2/3 for the odd and even states, respec- 
tively. Only positive values of the above expressions for 
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r 2 are of interest for the BI violation, since for t z < 
det(T) > 0, and the BI violation is impossible (see Sec. 

lira . 

For the both states we obtain [54[ , in view of Eq. ([(JTj) , 



S+ = 2 max <^ V2A a A fc sin 2/3, J \% A 2 sin 2 2/3 + t 



(98) 

Correspondingly, the optimal configurations for the odd 
and even states (with t z > 0) are given in Sees. IMG II 
and lIIIG~2| respectively. (The optimal configurations for 
the case t z < can be obtained, if necessary, from the 
results of Sec. IIIII ) In particular, the optimal configu- 
rations lie in the horizontal (a vertical) plane if the first 
(second) term in the braces in Eq. (|98|) is greater than 
the other term. 

According to Eq. (|98|) . for a given odd or even state 
the value of S+ depends on (3. In particular, when the 
first term in the braces in Eq. (f9"8"]l exceeds the second 
one, S+ = 2\/2A a Afc sin 2/3. This holds only when both 
A a Ab and /3 are sufficiently large, which implies relatively 
weak pure dephasing [14J . We do not dwell here upon 
the case of an arbitrary /3, since we are interested in a 
state which maximizes S+. 

In the ideal case, S max is obtained for a maximally 
entangled state [42|], but in the presence of decoherence 
this is not necessarily so. It is still of interest to consider 
S+ for the maximally entangled odd and even states (/3 = 
7r/4), i.e., \ 1 $>+) and |$+), respectively. For |^ + ) and 
|$+) Eq. (|9"8j) with (3 = 7r/4 becomes, respectively, 



S+ = Sl = 2[A 2 A 2 + max{A 2 A 2 , ( 7+ - l) 2 } 



1/2 



S + =Sl= 2[A 2 A 2 + max{A 2 A 2 , (1 - d) 2 } 



1/2 



(99) 



(100) 



When the first term in the braces in Eqs. (|9"9"|) and (|100|) 
is greater than the second term, Eqs. (|99|) and (|100[) 
reduce to Eq. (|93|) . One can show that 1 — d > | 7 + — 1|, 
the equality being for 7a = 7 (, equal to 1 or 0. This 
yields that Sf > S 1 * when the second term in the braces 
in Eqs. (|99[) and (|100[) is greater than the first term 
and t > 0. Hence, in the presence of decoherence, the 
BI violation for |$+) is generally greater than for |\I/+). 
This means that the BI violation duration re, defined by 
S+(t = tb) = 2, is generally greater for |$ + ) than for 

Generally, maximally entangled states are not optimal 
for the BI violations, unless the case (j93|) is optimal. Con- 
sider the maximal BI violations for the classes of the even 
and odd states. To this end, S+ in Eq. (|98p should be 
maximized with respect to (3. This yields for the odd 
state dnj 

S + = 2A a A 6 [l + max{l, ( 7+ - 1) 2 /(A 2 A 2 - 1 2 _)}] 1 ' 2 



if A 2 A 2 , >d u 
S+ = 2max{v2A a Af,,(i2} 



ifA?A?<di, 



(101a) 
(101b) 



where d\ = 7 2 + |t— |(t+ — 1) < 1 and di 
1| < 1, and for the even state (|75|) 



I7+ + I7-I 



•S-4 



2 A„ Ah 



1 + max < 1 



{1-df 



A 2 A 2 



d 2 



1/2 



if A 2 A 2 



> d, 



5+ = 2 



if A 2 A 2 < d. 



(102a) 
(102b) 



Equations ||101a|) and l|102a|) . in contrast to Eqs. (|101b[) 
and (|102b|) . can describe a violation of the BI. When the 
first term in the braces in Eq. (|101a|) or (|102a|) is greater 
than the second term, the case (IM1) is realized and the 
corresponding state is maximally entangled (/3 = 7r/4). 
In the opposite case, the optimal odd and even states 
generally are not maximally entangled, being character- 
ized by the following values of /3, respectively, 

/3 = arccos[ 7 -(7+ - 1)/(A 2 A 2 - 7 2 )]/2, (103) 
(3 = tt/2 - arccos[(d - d 2 )/(A 2 A 2 - d 2 )]/2. (104) 



Note that Eq. (|101b[) . with the second term greater than 
the first term, and Eq. (|102b[) are obtained for a nonen- 
tangled initial state: 1 10) if Tf < j| o r 1 10) if Tf > rf 
for Eq. (flOTbl) and |00) for Eq. pME|) . 

Numerical calculations show that the maximal BI vio- 
lation in the optimal even state is greater than or equal 
to that obtained in the optimal odd state. The differ- 
ences in S+ for the even and odd states can appear only 
in the case when the detectors are in a vertical plane. In 
this case the optimal states are generally nonmaximally 
entangled. The reason for this is that the observables 
whose axes do not lie in the horizontal plane are sensi- 
tive to energy relaxation. As a result, for instance, in 
the optimal even state the amplitude of 1 11) is less than 
the amplitude of |00) , since this bias reduces spontaneous 
decay and hence diminishes the detrimental effect of re- 
laxation on the BI violation. By the same reason, in 
the optimal odd state the amplitude of the excited qubit 
with a smaller 7 ^ (shorter Tf ) is reduced. However, in 
the case 7a = 7 & the optimal odd state, in contrast to the 
even state, is maximally entangled, and no relaxation re- 
duction is achieved. 

The states (fTTj) and (|T3"j) discussed here are special 
cases of the general odd and even states, respectively 
(0 < a < 2tt, < f3 < tt/2), 



e M sin/3|01), 



I*) =cos/3|10) 

|$) =cos/5|ll)+e M sin/3|00) 



(105) 
(106) 



The state (|105l) is relevant for experiments with super- 
conducting phase qubits [34], [36| . The states (|105[) and 
(I106[) result from Eqs. (|7T|) and (|73|) . respectively, under 
the rotation of qubit a around the z axis by the angle a. 
Hence, the results for 5*+ obtained in the present paper 
for the odd and even states hold, respectively, also for 
the states f|105[) and (|106p with the same (3, the corre- 
sponding optimal configurations being modified by the 
rotation of the detectors for qubit a around the z axis by 
a [see also the paragraph containing Eq. 
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E. Numerical results and discussion 

1. Pure dephasing 

First, let us discuss the case of the absence of energy 
relaxation, j a = i — 1, when decoherence occurs due 
to pure dephasing. Our calculations show that now |^+) 
and |3>+) are optimal states providing S max . In this case 
7+ = 2, 7_ = d = and both Eqs. (fTUTjl and (fT02]) yield 

GEES 



& 



= 2 



1 + A 2 A 2 . 



(107) 



Thus, now the BI can be violated for any level of deco- 
herence. The value (|107[) is achieved for the observation 
axes lying in a vertical plane [16j |. They are described 
by the results of Sees. IMG II and IIII G 2| taking into 
account that now T = diag(A a Ab, ±A a Afc, Tl), the upper 
(lower) sign correspomding to the odd (even) state. Note 
that the states (|105[) and (|106|) are also optimal. How- 
ever, not all maximally entangled states are optimal now 
[17| . since S is generally non- invariant with respect to 
local transformations. 



2. Identical decoherence of the qubits 

Next, consider the case of identical decoherence for the 
qubits of a pair. Now (i a = (lb = (i and 

la = lb = 7, A a = A h = A, d = 27(1 - 7), 7+ = 27. 

The numerical calculations show that in this case the 
maximal S can be always obtained with the even state 
(|75)l . i.e., iS ma x is given by Eq. (|102p with the account of 
Eq. (Tim 



S'n 



2X 2 



1 + max < 1 



On 



A 4 -d 2 



2 if A 4 < d. 



1/2 



if A 4 > d, 



(109a) 
(109b) 



Now the optimal odd state is the Bell state | V I / +), so that 



Eq. ([TUT]) reduces to Eq. 

trim 



(j9"9"j) with the account of Eq. 



S+ — S, 



2VA 3 



+ max{A 4 ,(2 7 -l) 2 }. 



(110) 



Figure [5] shows the dependence of 5+ on 7 with (i = 1 
(no pure dephasing) and /x = 0.9 for the even state 
(when S + = 5 max ) and the odd state |^+). Note that 
the straight segments in Fig. [2] (in particular, the both 
dashed lines) correspond to S+ = S h = 2\/2A 2 [Eq. ((9^)) ]. 
Equation (|109a[> and Fig. [2] imply that the violation of 
the BI is possible only for 7 > 2/3, this limit being ap- 
proached in the absence of pure dephasing (/1 = 1). As a 
result, for a given value of T±, the maximal BI violation 
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FIG. 2: S+ versus the relaxation parameter 7 for /1 = 1 (no 
pure dephasing) and /1 = 0.9. Solid lines: S+ — S ma x (|109[) 
with the even state (|73|l . dashed lines: Eq. (|110|l for the odd 
state |*+) (fT2|) . 



duration r^ is obtained for the even state when fi = 1, 
being given by 7 = e~ TB > Tl = 2/3 or tb = Tilnl.5 ~ 
0.405Ti. For comparison, we mention that in the case 
(i=l the odd state yields S + = Sh = 2-y/2~7, so that the 
BI can be violated only for 7 > l/\/2 w 0.707. This cor- 
responds to the longest tb for the odd state with a given 
Ti equal to [13, El tr = T\ In 2/2 w 0.347Ti. Note that, 
in contrast to Ref. |T(J, we obtain different values of tb 
for the even and odd states. 
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FIG. 3: S+ versus the pure dephasing parameter /x for 7 = 1 
(no decay) and 7 = 0.9. Solid lines: S+ = Smax (1109|) with 
the even state (|73|) . dashed lines: Eq. (|110|1 for the odd state 
l*+) Gill- 



Figure shows the dependence of S + on the pure de- 
phasing parameter (i. For 7 = 1 the plots for the odd 
and even state coincide and are given by Eq. (|107p . 
S + = Smax = yl + (i 4 = \1 + A 4 . In this case viola- 
tions of the BI can be achieved for any degree of pure 
dephasing. 

Figure|5]is a contour plot of S+ as a function of 7 and (i. 
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FIG. 4: Contour plot of the maximum S+ of S versus 7 and 
/i. Solid lines: the even state ((73)) (S+ = S ma x), dashed lines: 
the odd state |*+). 



The boundary of the region of the BI violation is shown 



by the solid line with 5 + = S n 



2. This boundary is 



obtained when the second term dominates in the braces 
in Eq. I|109a[) . The kinks on the curves in Figs. 121141 
correspond to a change of the dominating term in the 
braces in Eqs. I|109a|l and (fTT0| . Figures HH show that 
S+ for the odd state is generally lower than S max , which 
results in more stringent conditions on the decoherence 
parameters required for the BI violation than for the even 
state. The difference is significant when S + — 2 is small. 
However, for S+ > 2.4 there is practically no difference 
in the values of S+ for the odd and even states. 



3. No decoherence in one qubit 

Consider the extreme case of nonequal decoherence of 
the qubits, i.e., the case when decoherence is absent in 
one of the qubits, e.g., in qubit b. Now 75 = A& = 1, 
7+ = 1 + l a ; d = di = -7_ = 1 - 7 a , d 2 = 1, and 
hence Eqs. (jlOip and (|102[) coincide, i.e., both odd and 
even states give the same maximal BI violation S+. This 
S+ is maximally possible, S+ = S ma x, as shown by our 
numerical calculations, so that 



S'max = 2A Wmax^ 2, 



If A a > 1 - 7a, 



A Q -(l-7a) 2 



S+ = 2 if Xl < 1 - 7a . 



(Ilia) 
(111b) 



Now the optimal value of (3 follows from Eq. (|103[) or 
(HDD) to be 

= tt/2 - arccos{( 7a - 7')/^ - (1 - la) 2 }}/2. (112) 

Since qubit b is not affected by decoherence, S+ is in- 
variant with respect to arbitrary rotations of qubit b. 
Therefore, all states obtained from the optimal odd (or 



even) state by rotations of qubit b are optimal. [This ex- 
plains why the states (fTTj) and (|75|) yield the same results: 
these states are related by the unitary transformation a h x 
of the qubit b.] Moreover, since maximally entangled 
states transform to each other by a rotation of one qubit 
[44J (see Sec. MI Fi b). S+ is the same for all maximally 
entangled states [cf. Eq. (SS]| or (flUO"]) ], 



S + = 



/ 2V2A a , Tf < T 2 Q , 



x^vw 



2 T a > T a 
la ' J 1 — J 2 ■ 



(113) 
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FIG. 5: Contour plot of the maximum S+ of S versus 7a and 
fi a - Solid lines: S+ = S'max (|111[) for the even state (|73p and 
the odd state (|71[) . dashed lines: Eq. (|113[) f° r a maximally 
entangled state. 

Figure [5] shows the contour plot of S+ versus 7 a and 
[i a for the even and odd states which produce S+ = S'max 
(the solid lines) and for any maximally entangled state 
(the dashed lines). Equation (|111[) and Fig. [5] imply 
that now the BI violation is possible for 7 Q > 0.5, this 
limit being approached in the absence of proper dephas- 
ing (/i a = 1). Hence, in particular, for fi a — 1 the BI 
violation duration maximized over all states is given by 
t b = Ti In 2 w 0.693Ti. Equations (TlTTaj) and (TTTgj) im- 
ply that the BI violation for maximally entangled states 
is maximal (S+ = S'max) when there is no energy re- 
laxation (7a = 1) or when pure dephasing is sufficiently 
weak, so that the the first term in the braces in Eq. Q113P 
is dominating (see Fig. [5]). Note that maximally entan- 
gled states produce practically the same BI violation as 
the optimal states for 5+ > 2.4. 



4- General case 

In the general case when the decoherence parameters 
for the two qubits are different, the maximum of the BI 
violation can be obtained in a state of the form (|55|) . as 
discussed in Sec. IIVCI We performed several hundred 
calculations of S + for the states (|85)l and (|73]) with ran- 
dom values of the four parameters 7^ and fj,k (k = a,b) 
from the interval [0.8, 1]. For this parameter range the 
values of S ma x are greater than 2. In our calculations 
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•Smax resulted from the even state (l73|) in about 70% of 
the cases. In the cases, where the even state did not 
yield the maximal S, there were various optimal states, 
which included both general maximally entangled states 
(|14p and nonmaximally entangled states. However, the 
difference between S max and S+ due to the optimal even 
state was less than 0.1. 

Thus, there can be several approaches of a various de- 
gree of complexity and accuracy in order to obtain 5 max 
and the optimal observation conditions for given deco- 
herence parameters: (a) One can use the exact numerical 
approach of Sec. IIVCI which provides S'max, the optimal 
state, and the optimal detector configurations, (b) A sim- 
pler approach is to consider only the optimal even state 
[Eqs. (|102|) and (|104j) ]. which provides rather accurate, if 
not exact, result, as discussed above, (c) The analytical 
formula (|101|) and (|103|) for the optimal odd state can 
be used, if, e.g., in the experiment the odd state is re- 
alized more conveniently than other entangled states, as 
is the case for experiments with superconducting phase 
qubits. (d) An even simpler approach is to use the Bell 
state |* + ) or |$ + ) [see Eqs. flM} and flOO])]. In the ap- 
proaches (c) and (d) the BI violation is obtained under 
more restrictive conditions than in (a) and (b). However, 
if one requires a significant degree of the BI violation, say, 
S+ > 2.4 (which may be needed in the presence of other 
experimental errors), there is practically no advantages 
for one approach over the others (see Figs. [4] and [5]). 



DECOHERENCE AND MEASUREMENT 
ERRORS 



{Fim\ obeys F fc TO + F^ m = 1 and hence has two inde- 
pendent parameters, e.g., F$ — Fq and Ff = Ff 1; the 
measurement fidelities for the states |0) and |1), respec- 
tively. Equation (|114[) implies that in the presence of 
local measurement errors 

S = Tr(Bp), (115) 

where the error modified Bell operator [18| 

B = AB-AB' + A'B + A'B'. (116) 



Here, e.g., A = Q^(a) - Qg(S) [cf. Eq. ©] or A = 
C + £+a • ff a and B = £ b _ + Qb ■ a b , where 

& = F$ + F?-l, g_=Ff-F$ t (117) 



A' and B' following from A and B on the replacement of 
a and b by a' and &', respectively. 



B. Modified Bell operator for decoherence and 
errors 

In the presence of decoherence and errors we perform 
the substitution (|74[) in Eq. (|115p . It is useful to recast 
the resulting expression as 



S = Tr(Sp ) 



(118) 



where the Bell operator modified by errors and decoher- 
ence is 



In the previous sections we assumed that measure- 
ments are ideal. Here we take into account the possibility 
that measurements of the qubits are performed with local 
(independent) errors. Effects of local errors were studied 
elsewhere [18l |55[. In this section we discuss combined 
effects of local errors and local decoherence. The present 
case is rather involved, since it includes complications 
due to both decoherence and errors. Here we discuss 
only the general approach to the problem, whereas a de- 
tailed analysis is out of the scope of the present paper 
(see also Ref. UM . 



A. Description of errors 

In the frame of the error model considered in Ref. [l| 
the measured probabilities are written in the form 



Pff 



E 



^im^jnPmn — ^(QiQjP): 



(114) 



where Q\ = YL=Q F L,Pm (* = 0,1; k = a,b) arc the 
POVM operators describing the measurement Q [cf. Eq. 
0] and F* m is the probability to find qubit k in state 
\i) when it is actually in state \m). The fidelity matrix 



B = {CI ® C* b )(B) =AB- AB' + A'B + A'B 1 . (119) 

Here A = Qt{a) - Qg(S), B = Q\{b) - Q b Q (b), etc., where 
the POVM operators Q\ are given by Q\ = C* k (Qi) (k = 
a, b) or [cf. Eq. (S3])] 

j — 1 771=0 j = l 

(120) 
so that [cf. Eq. 



A = e+C+ha-l)a z +e+qa-Va, (121) 

whereas the operators A', B, and B' are given by Eq. 
(I12ip . where a is replaced by a', b, and b', respectively. 

Consider properties of S which can be helpful in cal- 
culations. One of the peculiarities present in the case 
with errors is that the maximal and minimal values of S 
for a given state are generally not equal by the magni- 
tude [13| , S+ 7^ 1-5-1. Note, however, the relations which 
follow from Eqs. (118). (119), and (121), 

S -> -S if a -v -a, a' -v -a', F a *-» Ff; (122a) 
S^-S if b -> -b, b' -» -b', F b *-> F b . (122b) 
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Moreover, one has S+ = \S-\ if the two measurement 
fidelities are equal, at least, for one qubit: 



F a = F? or F b = F?. 



(123) 



The relations (|122p and (|123[) were obtained previously 
(18j for the case without decoherence. 

In the special case of equal measurement fidelities for 
each qubit, 

F a = Ff = F a , F 6 = F* = F b , (124) 

Eq. pT9]) yields B = (2F a - l)(2F b - l)B d and hence 

S = (2F a -l)(2F b -l)S d , (125) 

where Sd is the value of S obtained in the presence of 
decoherence but in the absence of measurement errors. 
This case can be analyzed, as discussed in Sec. IIVI 

C. Discussion 

There is no analytical solution in the presence of er- 
rors [lg, [55|. Moreover, the present case involves an 
eight-dimensional parameter space (there are two deco- 
herence parameters and two measurement fidelities for 
each qubit), which additionally complicates the analysis. 

Similarly to Sec. IIVI the modified Bell operator B 
(I119P can be used for numerical calculations, since for 
given decoherence and error parameters the maximum 
(minimum) value of S equals the maximum (minimum) 
of the greatest (smallest) eigenvalue of B over the detec- 
tor directions, the optimal state being given by the corre- 
sponding eigenvector. Now, as in Sec. IIVI S is invariant 
only to rotations of the qubits and detectors around the z 
axis, which allows one to reduce the number of the fitting 
parameters from eight to six by setting, say, a y = b y = 0. 
This number cannot be further reduced, since, in contrast 
to Sec. IIV| in the presence of errors the set of optimal 
configurations for a given state is not continuous 18], 
and hence it is impossible to choose one of the optimal 
measurement axes at will. 

This computation procedure is relatively slow. It pro- 
duces generally different optimal states for different val- 
ues of the parameters. An approach which is faster and 
more relevant for most experiments is to consider the 
BI violation for a specific initial state po, e.g., the odd or 
even state [Eqs. (|7T|) and ([73]) ]. In this case an expression 
for S resulting from Eq. (|118[) is varied over the detector 
directions and perhaps the state parameters [e.g., (3 in 
Eqs. CLO) and ([73])]. 

Note that in the cases of the odd and even states the 
number of the detector parameters can be reduced from 
eight to seven. Indeed, since the odd (even) state is in- 
variant under a rotation of qubits a and b around the z 
axis by an arbitrary angle a (a and —a, respectively), 
one can set, e.g., a y — 0, in view of the z-rotation sym- 
metry for S. For a detailed analysis of the case of the 
odd state, Eq. (J7TJ) , see Ref. GH 



VI. CONCLUSIONS 

In the present paper we have considered conditions for 
maximal violations of the Bell inequality in the presence 
of decoherence. In addition, combined effects of decoher- 
ence and local measurement errors have been discussed. 

Since decoherence transforms a pure entangled state 
into a mixed state, we have begun the consideration from 
the study of optimal conditions for the violation of the BI 
(TIJ for a general (pure or mixed) state. We have obtained 
all detector configurations providing the maximal value 
of the parameter S in Eq. (fT]) for an arbitrary state. We 
have shown that generally the set of all optimal configura- 
tions for a given state is characterized by one continuous 
and one discrete parameters, whereas in special cases it 
can be characterized by two or three continuous param- 
eters. We have obtained also the symmetry relation for 
the optimal detector orientations, Eq. ([50]) or (|5"T)) . which 
follows from the invariance of S with respect to the qubit 
swap. 

Further, we have considered effects of local decoher- 
ence on the BI violation. The decoherence model used 
includes energy relaxation at the zero temperature and 
arbitrary dephasing. We have employed the decoherence 
superoperator in the operator-sum form with the Kraus 
operators which generalize the previously known ones. 
We have expressed S as the average over the initial state 
of the Bell operator modified by decoherence. This op- 
erator has been used for numerical calculations in order 
to obtain the maximal BI violation for any values of the 
decoherence parameters. We have also developed a faster 
numerical approach, which is based on an analytical so- 
lution for a certain class of states. 

We have studied the BI violation maximized over the 
observation directions and over all states or a class of 
states and obtained all corresponding optimal detector 
configurations. We have obtained simple analytical solu- 
tions for the odd and even states [Eqs. (|7T|) and ([75]) or, 
more generally, (|105[) and (|106[) ]. which are often used 
in experiments on the BI violation. In particular, the 
odd state is relevant for experiments with superconduct- 
ing phase qubits. Whereas in the absence of decoherence 
the optimal detector configurations for the odd and even 
states are vertical (Sec. IIIIG3|) . in the presence of deco- 
herence the optimal configurations are either vertical or 
horizontal. We have discussed both the general case of 
arbitrary decoherence parameters and a number of im- 
portant special cases. In particular, we have discovered 
that the even state is optimal in most cases. Our analysis 
have been illustrated by numerical calculations. 

We have considered the combined effects of local errors 
and decoherence. In this case the maximal Bell violation 
depends on eight parameters. We have derived the Bell 
operator modified by decoherence and errors and have 
used it to discuss symmetry properties of S and special 
cases. We have also outlined the numerical approaches 
in this case. 

We have discussed the application of the above analysis 
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to superconducting phase qubits. However, the present 
results are applicable to many types of qubits. Moreover, 
the present results have relevance to the ongoing discus- 
sion of effects of decoherence on entanglement, a major 
resource in the field of quantum information. 
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